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Abstract. Geometric optics is analysed using the techniques of Presymplectic 
Geometry. We obtain the symplectic structure of the space of light rays in a medium 
of a non constant refractive index by reduction from a presymplectic structure, and 
using adapted coordinates, we find Darboux coordinates. The theory is illustrated 
with some examples and we point out some simple physical applications 
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1. Introduction 



The Hamiltonian formulation of Classical Mechanics is introduced in most 
textbooks as a method of transforming the set of second order Euler-Lagrange 
differential equations into a first order system (see e.g. [1]). However, systems 
described by Hamilton-like equations are interesting not only in Physics but also 
in Mathematics, Chemistry, Biology, etc., and this caused a big interest in the 
study of these systems. Moreover, the introduction in a geometric way of the con- 
cept of Hamiltonian Dynamical System [2], which reduces in simple cases to the 
Hamiltonian or the Lagrangian formulation of the Mechanics, has motivated the 
development and study of symplectic structures and related tools. This enables us 
to deal, in a coordinate free way, with mechanical systems for which the configu- 
ration space is not topologically trivial or even it is an infinite-dimensional linear 
space as it happens in the case of Field Theory. 

The Hamiltonian treatment of Geometric Optics, based on Fermat's princi- 
ple, preceded its applications to Classical Mechanics, based on the corresponding 
Hamilton principle. This suggested that a Hamiltonian formulation of Geometric 
Optics using the tools of Modern Differential Geometry may also be useful [3]. 
Then, this geometric formulation of Optics is receiving much attention during the 
last ten years and it is now well established. For instance, a recent book on Math- 
ematical Methods in Physics [4] includes a Chapter IX on matrix and Hamilton 
methods in Gaussian optics. In particular, the development of Lie methods in 
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Optics [5-10] has been very fruitful and allowed a new perspective in the theory 
of aberrations, where the symmetry theory imposes "selection rules" for the aber- 
ration coefficients [11-13]. The main difference of Geometric Optics with particle 
mechanics is the existence of refracting surfaces separating regions of a different 
constant refractive index. On the other hand, there is no natural "time parame- 
ter" for the ray and it is not the parametrized curve but the trajectory itself, no 
matter its parametrization, what is actually relevant in Geometric Optics. 

The usual geometric approach to Symplectic Optics is based on the fact that 
an optical system can be seen as a black-box relating constant refractive index re- 
gions and then incident and outgoing rays are characterized by appropriate points 
(<ZiiuPin) an d (<?out,Pout) of a phase space in such a way that the optical device can 
be considered as a canonical transformation in such phase space. This reminds 
us what happens in Quantum scattering theory where |in) and |out) states are 
related by a unitary transformation. 

It is very often said in physics textbooks that active and passive viewpoints 
in transformation theory of physical systems are equivalent. However, this is 
only true for the simplest case in which the manifold describing the system is M. n . 
Actually, coordinates can only be introduced in a local way while a transformation, 
from the active viewpoint, is a global concept. So, if M is a differentiable manifold, 
the configuration space of a system, the phase space, or any other similar thing 
else, a transformation of M is a diffeomorphism F : M — > M. It induces a change 
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of charts, the chart (U, 4>) becoming (W, <fi') with W = F(U) and <fi' = (poF -1 . On 
the contrary, a change of coordinates is a local concept and it does not produce 
any diffeomorphism of the manifold M. 

We will first prove that actually the space of light rays in a constant index 
medium is a symplectic manifold and that Darboux coordinates for these rays are 
the q and p usually chosen. Moreover, it will be shown using the techniques of 
the geometric approach to singular Lagrangian systems that, under very general 
assumptions, the space of light rays in an arbitrary medium is a symplectic man- 
ifold, and if we restrict ourselves to optical systems such that the refractive index 
depends only on the third coordinate x 3 and the index takes, possibly different, 
constant values for x 3 > L and x 3 < L, we can choose Darboux coordinates by 
fixing the x 3 coordinate in any of these two regions and taking Darboux coordi- 
nates for the corresponding problems of constant index. This property provides 
a justification of the choice of coordinates as usually done for the ingoing and 
outgoing light rays in the corresponding constant index media. Therefore, it is the 
choice of "Darboux coordinates" in the ingoing and outgoing regions where light 
moves freely (i.e., in a constant refractive index medium) what motivates that the 
passage of light through the optical system seems to be a canonical transformation 
(active point of view) when actually it only corresponds to a change from one to 
another Darboux coordinate system (passive viewpoint) in a symplectic manifold. 

The paper is organized as follows: For the sake of completeness we give in 



Symplectic structures in Optics 5 
Section 2 a short summary of important concepts of modern Differential Geometry, 
which are well konwn only for few opticists, and introduce the notation to be used 
in the paper. In Section 3 we develop an idea introduced in [14] to show in a 
geometric way that the space of oriented straightlines in the Euclidean three- 
dimensional space M 3 is the tangent bundle TS 2 of the two-dimensional sphere, 
and then the identification with T*S 2 , given by the Euclidean metric, shows that 
this space is endowed in a natural way with a symplectic structure. We will find out 
Darboux coordinates, which essentially reduce to those used by Dragt, Wolf and 
coworkers for the description of light rays in a constant refractive index medium. 
In Section 4 we show how the minimal optical length Fermat's principle leads to a 
symplectic formulation of Geometric Optics, by making use of the singular optical 
Lagrangian and its relation with a regular Lagrangian. The reduction theory of 
presymplectic manifolds, following the ideas developed by Marsden and Weinstein 
[16], will be used, and we will illustrate the method for finding coordinates adapted 
to the distribution defined by the kernel Kerc^L of the presymplectic structure 
defined by the singular optical Lagrangian in the case of a constant index medium 
in both cartesian and cylindrical coordinates. The most general case of a variable 
index medium is considered in Section 5 and the complete solution is obtained 
for the particular cases of systems in which either the index n depends on the 
third coordinate x 3 alone or the very interesting case in which the system is axis- 
symmetric and n is a function of the distance to this axis [17]. Once Darboux 
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coordinates have been found we can consider the problem from the active viewpoint 
and take advantage of the algebraic methods recently developed for computing 
aberrations. 

2. Notation and basic definitions 

The existence of constraints for mechanical systems motivates that the config- 
uration space is not anymore but some subset Q of points, called configuration 
space. Physicists know the convenience of using generalized coordinates for the 
description of such systems. The corresponding geometric concept is that of dif- 
ferentiable manifold. We recall that a chart for a topological space M is a pair 
(U, (f>) where U is an open set of M, <j>(U) is an open set of R n and : U — > (f>(U) 
is an homeomorphism. Two charts (U, <fi) and (U', </>') are said to be compatible 
if either U n W = or cp' o : <f>{U n W) -> <f>'(U n W) is a diffeomorphism 
of open sets in IR n . A differentiable structure in M is given by an atlas, a set 
of compatible charts covering M. From the intuitive point of view that means 
that there is a way of parametrizing points by domains in M and coordinates for 
overlapping domains are related by differentiable expressions. The set M endowed 
with a differentiable structure is said to be a differentiable manifold. We will say 
that a function / : M — > R is differentiable in a point m E M when there is a 
chart (U,4>) such that m G U and / o is differentiable. In a similar way a 
function F : M — > AT between differentiable manifolds is differentiable at a point 
m 6 M if there exist charts (W, 0) in m and (V, ip) in F(m) G such that the 
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map ipoFocj) -1 is differentiable. An invertible differentiable map F such that F _1 
is also differentiable is called diffeomorphism and a diffeomorphism F : M — > M 
corresponds to the physical idea of a transformation of M, in the so called active 
point of view. However such a diffeomorphism produces a changes of coordinates, 
the chart (U t <f>) becoming (U',<f>') with W = F(U) and <p' = <p o F -1 . This is 
called the passive viewpoint. An expression like y l = F l {x^) can be seen either as 
a transformation in which the point of coordinates y % is the image under F of the 
point of coordinates x\ or alternatively as a change of coordinates, from x % to y l , 
for the same point. 

Vectors in a point m G M may be introduced by an equivalence relation among 
curves starting from m. So, if 71 (0) = 72(0) = m are two such curves they are 
equivalent if 071 is tangent to ^072 at 4>{m). They can alternatively be seen as a 
map assigning to each differentiable function in m the real number o 7i)| t= o, 
and then having in mind that for any curve 7, 

|*=0 V UX J \4>{m) 



the corresponding vector will be denoted v = v 1 ^^ with v l = aKq> 2t" \t=o~ Then, 
if 4> o 7! is tangent to o 72, f t {f o -yi)|t=o = f t U 72)|t=o- The set of a11 vectors 
in the point m G M is a linear space called the tangent space at m and denoted 

T m M. 

The set of all possible vectors in all points of M is called the tangent bundle 
TM. It is endowed with a differentiable structure obtained from that of M . In 
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fact, a basis of T m M is obtained by taking the vectors tangent to each coordinate 
line. The map r : TM — * M assigning to each vector the point where it is placed 
is such that for any chart in M, r _1 (W) = Wx R n . The corresponding coordinates 
will be denoted (q l ,v l ). A similar process can be followed by glueing together 
the dual spaces T^M and we will obtain in this way the so called cotangent 
bundle n : T*M — > M, with coordinates induced from base coordinates to be 
denoted (q l ,Pi)- It is remarkable that it is possible to use in both spaces TM 
and T*M more general ways of introducing coordinates, mixturing base and fibre 
coordinates. However the ones introduced previously are adapted to the tangent 
or cotangent structure, respectively, while the more general coordinates will hide 
that structure. Sections of these bundles are called vector fields and 1-forms in 
M, respectively. 

The important point is that tangent and cotangent structures are characterized 
by the existence of canonical objects. The best known one is a canonical 1- 
form 9 in T*M that in coordinates (q l ,Pi) adapted to the cotangent bundle looks 
6 = pi dq l . On the other side, the tangent bundle is characterized by a (1, l)-tensor 
field called vertical endomorphism S that in terms of natural coordinates (q\v l ) 
of the tangent bundle TM is given by 



The exterior differential of the canonical 1-form 9 in a cotangent bundle is a 
closed 2-form u = — d6 such that uj An ^ 0, and then this nondegeneracy allows 




<g) dq\ 



(2.1) 
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us to put in a one to one correspondence vectors and covectors in a point, and by 
extension, vector fields with 1-forms in M, in much the same way as it happens 
with Riemannian structures in General Relativity. In the natural coordinates, 
u = dq 1 A dpi. The generalization of this structure is the concept of symplectic 
manifold, a pair (AT, fi) where O is a nondegenerate closed 2-form in the differen- 
tiable manifold N. The dimension of iV is an even number 2n. Moreover, cotangent 
bundles are the local prototype for this more general objects, as established by 
the well known Darboux Theorem: given a point m G N in a symplectic manifold, 
there is a coordinate neighborhood of m with coordinates (y 1 ,... ,y 2n ) (called 
Darboux coordinates) such that Q is written Q = dy l A dy n+l . Summation from 1 
to n for the index i is understood. It is usual to denote these new coordinates as 
y l = q % and y n+l = Pi, for % = 1, . . . , n, for which the expression looks like for the 
symplectic T*M manifold. Of course these Darboux coordinates are not uniquely 
determined. For instance we can introduce new coordinates Pi = —q 1 , Q l = Pi, 
that are also Darboux coordinates. A change of base coordinates Q l = Q l (q) in 
a cotangent bundle induces a change in fibre coordinates pi = PjdQi /dq 1 in such 
a way that (Q\ Pi) are also Darboux coordinates. Notice that in the general case 
of an arbitrary symplectic manifold the notation of g's and p's is arbitrary while 
in the T*M case g's are actually coordinates on the base space and p's are fibre 
coordinates and therefore these last ones should take all possible real values. Now, 
the choice of an arbitrary function H E C QO (N) 1 determines the vector field Xh 
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such that i(X H )uj = dH, i.e., u(X H ,Y) = YH, VY G 3t(N), that in Darboux 
coordinates is 

X H = — — -— — (2.2) 
dpi dq i dq i dpi ' 

Its integral curves will be determined in Darboux coordinates by the well-known 
Hamilton equations. 

There is no natural symplectic structure in TM, but for any function L e 
C°°(TM), we may define an exact 2-form in TM, uj l = —dO^, with the 1-form 9 L 
being defined by Ol — dLoS, and a function E L = A(L)— L, called energy function. 
In the above mentioned coordinates of TM we have the following expressions: 

«L = ^dq\ (2.2) 

^ = ^^ Ad «' + &'«' A *"' (2 ' 3) 

A = "'^ < 2 ' 4 > 
E L = ^_-L, (2.5) 

Here A e X(TM) denotes the Liouville vector field generating dilations along the 
fibres. When is nondegenerate it defines a symplectic structure on TM, and 
a vector field Xl determined by i(Xl)u>l = dE^. The Legendre transformation 
TL : TM — > T*M relates the symplectic structures in TM and T*M respectively, 
and if H e C°°(T*M) is defined by if o TL = E L , then TL*X L = X H . 

Presymplectic structures may arise either when using some constants of motion 
for reducing the phase space or also when the Lagrangian that has been chosen is 
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singular. Then we will have a pair (A^ ,O ) where fi is a closed but degenerate 
2-form. A consistent solution of the dynamical equation can only be found in 
some points, leading in this way to the final constraint submanifold N introduced 
by Dirac (see e.g. [15]). The pull back Q of the form Qq on this manifold will be 
assumed to be of constant rank. The recipe for dealing with these systems was 
given by Marsden and Weinstein [16]. First, in every point m G N, Kerfi m is a 
fc-dimensional linear space, so defining what is called a /c-dimensional distribu- 
tion. The important point is that closedness of O is enough to warrant that the 
distribution is integrable (and then it is called foliation): for any point m G N ', 
there is a /c-dimensional submanifold of N passing through m and such that the 
tangent space at any point m' of this surface coincides with Ker fi m / . Such integral 
/c-dimensional submanifolds give a foliation of N by disjoint leaves and in the case 
in which the quotient space N = N/ Ker O is a differentiable manifold, then it is 
possible to define a nondegenerate closed 2-form O in N such that 7r*0 = O. Here 
n : N — > iV is the natural projection. It suffices to define Cl(vi,V2) = fi(vi,V2), 
where v\ and are tangent vectors to iV projecting under tt onto v\ and 
respectively. The symplectic space (N, O) is said to be the reduced space. 

3. The symplectic structure of the light 
rays space in a constant index medium 

In a recent paper [14] it was shown that the set of oriented geodesies of a 
Riemannian manifold can be endowed with a symplectic structure. We will study 
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here some particular examples of interest in Symplectic Optics. To start with, 
let us see that the set of oriented straightlines in the plane, that is well known 
to be the set of light rays in a two-dimensional constant rank medium, can be 
endowed with a symplectic structure. Moreover, such a set can be considered as 
the cotangent bundle of the one-dimensional sphere S 1 , i.e. a circumference. If 
an origin O has been chosen in the plane, every oriented straightline that does 
not pass through the point O is characterized by a unit vector s pointing in the 
line direction and a vector v orthogonal to s with end on the line and origin in O. 




Fig.l 



So, the straightlines of a pencil of oriented parallel lines have the same direction 
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and only differ in the distance to the origin, and therefore they are characterized 
by proportional vectors v and the same s. Straightlines passing through O with 
direction given by s correspond to v = 0. The vectors v and s being orthogonal 
and s.s = 1, the couple (s, v) can be seen as a tangent vector to the unit circle S 1 
at the point described by s as indicated in the following figure: 
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Fig.2 



The circle 5" 1 is a Riemannian one-dimensional submanifold of the Euclidean 
two-dimensional space and the Riemannian metric in S 1 can be used to identify 
in each point s the tangent space TgS 1 with its dual space T^S 1 and therefore the 
tangent bundle TS 1 with the cotangent bundle T*5' 1 . This identification shows us 
that the space of oriented straightlines in the Euclidean two-dimensional space can 
be endowed with an exact symplectic structure which corresponds to the canonical 
structure for the cotangent bundle T*S X . 

Let us look for the kind of Darboux coordinates for such a symplectic form 
mentioned in the preceding section. They will be induced from a choice of coordi- 
nates in the base space. As it is a circle, a good choice will be an angle coordinate. 
A straightline y = mx + b with slope m = tan 9 will be represented by a vector 
orthogonal to the vector s = (cos 6, sin 9), and length b cos 9, namely, 



On the other hand, the vector d/d9 is unitary in the Euclidean metric, and 
therefore the point (9,po) E T*S 1 corresponding to (9,vq) is given by pg = vg. 
Therefore, the symplectic form in TS 1 translated from the canonical symplectic 
structure in T^S 1 uj = d9 A dpg will be 




(3.1) 



lo = d9 A d(bcos9) = <i(sin6>) A db, 



(3.2) 
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which shows that Darboux coordinates for uj adapted to the cotangent structure 
are not only (9, bcos9) but also 

q = sin 9, p = b, (3-3) 

which are more appropriate from the experimental viewpoint. So, the flat screens 
arise here as a good choice for Darboux coordinates. 

The choice of the norm equal to one for the vector s representative of the line 
direction is arbitrary. In the case of Geometric Optics, the light rays in a medium of 
a constant index n are straightlines. For reasons which will become clear later, the 
choice usually done is s.s = n 2 , the Darboux coordinate q then being q = n sin6>. 
This leads to the image of the Descartes sphere, a sphere of radius n whose points 
describe the ray directions (see e.g. [20]). Therefore, the space of light rays in 
a medium of constant index is like a phase space of a system with configuration 
space the sphere S 1 of radius n corresponding to ray directions. It is noteworthy 
that it is very often used a somehow misleading notation of p for the product n sin 9 
even if the last one corresponds to a coordinate in the configuration space. In our 
identification it is clear that the coordinate q is bounded while p is unbounded, as 
it was expected to be for a cotangent bundle structure. 

The study of oriented straightlines in Euclidean three-dimensional space fol- 
lows a similar pattern, but the search for Darboux coordinates is a bit more 
involved. A straightline r will be characterized by a vector s, running a two- 
dimensional sphere, which fixes the ray direction, and a vector v orthogonal to r 
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with origin in O and end in r. The conditions s.s = 1 and s.v = mean that the 
couple (s, v) determines an element of the tangent bundle TS 2 . The points (s, v) 
of the zero section correspond to straightlines passing through the origin. 

The Riemann structure in S 2 coming from the Euclidean structure in R 3 pro- 
vides us the identification of TS 2 with T*S 2 and then it allows to endow in this 
way TS 2 , and therefore the set of oriented straightlines in M 3 , with a symplectic 
structure coming from the canonical structure in T*S 2 . In order to find Darboux 
coordinates for it let us consider, for instance, the local chart in S 2 with domain 
in the upper hemisphere and the map <p(x, y, z) = (x, y). In other words, the point 
in S 2 with local coordinates (1*1,142) is (111,112, a/1 — u 2 — itf). The corresponding 
basis for the tangent space in a point of coordinates (u±, 1x2) is 

/in IXi 1X2 x 

e Ul = (l,0, ), e U2 = (0,l, ) 

y/l-Uf- U\ y/l-u(- u\ 

and the Riemannian metric will be given in these coordinates by 

9(ui,u 2 ) = - 1 —( 1 - Ul ? 1U \). (3.4) 

v ; 1 - u\ - u\ \ uiu 2 l-u(J v ' 

Then, a vector v tangent to the sphere S 2 can be expressed either as v = 
v x e x + v y e y + v z e z or as v = a e Ul + b e U2 . 

In particular, it is easy to check that a = v x and b = v y . In the identification 
of T S S 2 with T*S 2 given by the Riemannian metric, the corresponding point will 
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have coordinates 

\ — u\ U\U 2 

\ — u\ — u\ 1 — u\ — u 2 

U\U 2 1 — u\ 

1 — uf — U2 1 — u\ — U2 
and therefore the symplectic 2-form in TS 2 will be given by 

, (1 - ul)v x + u 1 u 2 v y \ A ,(u 1 u 2 v x + (1 - u\)v 

u = dui A a ( ; 2 2 ) + d u 2 A a 1 



= du\ A (i 
+ ^2 A d 



1 — u( — u% J ' V \ — u\ — ui A 

(1 - npf s - u\v x + m| i> x + UiU 2 V y \ 
\-u\-u\ ) 

UlU 2 V x + (1 - )u y - u\v y + 

1 - u\ — u\ 



I Ui ~UlV x - U 2 Vy 

— du\ A a I v x 



+ du 2 A d \ v v — 



U 2 —U\V X — u 2 v 



y 



y/l-uj- u\ y/l-uf- 



2 



u 2 



•1 u l I / U 2 

= du\ Ad I v x . = =v z + du 2 A \ v v 



R 2 2 z y y R 2 



U 



2 



(3.5) 



a/1 - u\ - u\ J \ a/1 - uj - u\ 

that can be rewritten as u> = du\ A db x + du 2 A db y , with 



v * 7= U \ -> v z b y = v y u U \ > v *- ( 3 - 6 ) 



The geometric interpretation is clear. The point of cartesian coordinates 
(b x , b y ,0) is the intersection point of the straightline r with the plane it defined 

by z = 0. In fact the points of r are given by v + As for any real number AeM. 

v z 

The intersection point of r with tt corresponds to the value A = , and then its 

Sz 

coordinates turn out to be (b x , b y , 0). So, in the study of the set of light rays in a 
three-dimensional constant index medium, Darboux coordinates are, for instance, 
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given by the projections e x .s and e y .s, as (/-coordinates and the cartesian coordi- 
nates of the intersection point of r with n as momentum coordinates. Therefore, 
it also shows that flat screen are appropriate as a method of introducing Darboux 
coordinates for the symplectic form. 



Light rays trajectories in Geometric Optics are determined by Fermat's prin- 
ciple, according to which the ray path connecting two points is the one making 
stationary the optical length 



This corresponds to the well-known Hamilton's principle of Classical Mechanics 



with an "optical Lagrangian" L = n y v\ + v 2 + v%, which is a differentiable func- 
tion in TR S up to at the zero section. 

The point to be remarked here is that such Lagrangian is a homogeneous 
function of degree one in the velocities and consequently L is singular and the 
corresponding energy function vanishes identically. The singularity of L is related 
with time reparametrizacion invariance (see e.g. [21] for a geometric approach to 
the second Noether's Theorem). Moreover, it is possible to relate the solutions of 
the Euler-Lagrange equations for L with those of the regular Lagrangian L = ^L 2 , 
up to a reparametrization. The geometric explanation of this fact was given in [14] . 
We will study here the presymplectic structure defined by ul and the reduction 
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(4.1) 
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procedure giving rise to a symplectic structure in the quotient manifold. We will 
also find Darboux coordinates for the reduced symplectic form. For the sake 
of simplicity, let us first consider in this section the simplest case in which n is 
constant. 

The coordinate expressions of the 1-form 9l and the 2-form lvl are (see [14]) 
n 2 

9 L = — (v x dx + v y dy + v z dz) 
n 2 

ujl = —(dx A dv x + dy A dv y + dz A dv z ) (4-2) 

LJ 

n 4 

+ Jz( v * dv x + v y dv y + v z dv z ) A (v x dx + v y dy + v z dz), 
the kernel of ujl being generated by the Liouville vector field 

A = Vx ir +Vy ir +Vz ir (4 - 3) 

and the dynamical vector field giving the dynamics corresponding to L, 

d d d 

r = V *8- X +V »ey- + V *dz- (4 ' 4) 

The 2-form ujl is closed and defines an integrable two-dimensional distribution 
V generated by A and V. The relation [A, V] = T, expressing that the vector field V 
is homogeneous of degree one in velocities, shows explicitly that V is an involutive, 
and therefore integrable, distribution. 

In order to find coordinates adapted to V we can consider that, in a neighbour- 
hood of a point for which v z ^ 0, the distribution V is generated by the commuting 
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vector fields A and K = — V. The general theory of integrable distributions, see 
e.g [22], says us that there exist local coordinates (x±, X2, x$, yi, 2/2, Vi) such that 

A=A K= JL. (4 . 5) 
0x3 dy 3 

The coordinates x\,X2,yi and 1/2 will be given by functions solution of the 
system 

df df df 
A/ = v x h v y h v z — = 

dv x di> z 

1 / df df df\ 
Kf= - z { Vx d^ + Vy ^ +Vz d-z) =0 ' 



while X3 will be given by a solution of 



a p df df df 

A/ = v x h v y \~v z — 

OV x OVy dv z 



Kf = — [v x — + v y — + v z — =0, 

v z \ ox dy dz 



and ys by a solution of 



df df df 



A/ = V x h U„ h t> Z — = 

dv x dv y dv z 



, e I f df df df . , 

#Y = — ( ^- + u v -5- + -5- = 1- 

v z \ dx dy dz 



A particular solution of these systems is given by 

Vx V X Vy 

X\ = Z - X, X 2 = y X, x 3 = logv z 

v z v z v z 

V X Vy 

yi = — , V2 = — , V3 = z 

v z v z 
The inverse transformation of coordinates is 

X2 - V2Xl 

x = yiy3-xi, y = , z = y 3 

Vi 

v x =yie x \ Vy = y 2 e x \ v z = e X3 



(4.6) 



(4.7) 
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When written in these new coordinates the expression of the 2-form ujl be- 



comes 

ne 3x 3 

U) L = 



dyi A dxi H ^ 1 dx2 A eij/2 + —dy 2 A ctei 

2/i 2/1 



. x 2 -xiy 2 + x 2 yf 

H 2 "2/2 A dyi + y 2 ayi A cte 2 

2/i 

which can also be rewritten with a reordering of terms as 

UL = (i + j/ 2 + y 2 )3 / 2 i d (- x i) A ((! + 2/2)^2/1 - 2/12/2^2/2 



) 



2/1 



or 



^ W J I S/i \ , j / x 2 - 2/2^1 \ . , / 2/2 



v/i + 2/? + 2/1 y V 2/i / V V 1 + vl + vl 

(4.8) 



and when using the old coordinates, 



Vx , X2 - 2/2^1 , 

- xi = x z = b x , = y z = by 

v z 2/i v z 

nyi nv x 
y/1 + y\ + 2/ 2 2 ^J v l + v 2 + v 2 z X (4-9) 
n 2/2 _ n« 9 
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Py 



y/1 + y\ + y\ ^ V 2 + v 2 + 

from which we have got Darboux coordinates for the reduced form. Notice that the 
explicit expression of ujl in these coordinates shows that it passes to the quotient, 
because its coordinates do not depend on xs and j/3, and then there exists a 
symplectic form O in the quotient space with the same coordinate expresion as in 
(4.8) for ujl- The quotient space is but the space of light rays and therefore it is 
TS 2 as indicated before. 
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We can also use cylindrical coordinates (r, 6>, z) for solving the problem. The 

vector fields IY and A are then given by 

d d d 2 d n v r V0 d 
Tl = v r — + V0— + v z — + rv e - 2 



dr 86 dz dv r r dvg , A -, n \ 

a a a < 410 > 

A = V rx l-"9n ^ v *-^— 

OV r OV0 ov z 

and a similar computation, choosing now K' = ^-Tl and A as generators of V, 
leads to the following adapted coordinates 

I q V r l rV r v z , , 

£i = #-arctan , x 2 = z 5 5—5-, a; 3 = logv z 

rvg v£ + T 1 Vq 

1 2 v e , v 2 +r 2 vj , 
Vi = r—, y 2 = 2 — ' V3= r i 

and Darboux coordinates for the induced form in the quotient are 

v r rv r v z 
£1 = — arctan , £2 = z 



rvg ' v 2 + r 2 Vg 

nr 2 ve nv z 



(4.11) 



\/vf+ r 2 v 2 B + v 2 ' \Ar + r 2 i;^ + 

5. The case of a medium of a non-constant index 

Let us now consider the most general case in which the refractive index of the 
medium is not constant but it is given by a smooth function n(x 1 , x 2 , x 3 ). Fermat's 
principle suggests us to consider the corresponding mechanical problem described 
by a singular Lagrangian L{q,v) = [g(v, f)] 1//2 , where g is a metric conformal to 
the Euclidean metric go, 

g(y, w) = n 2 g (v,w). (5.1) 
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This problem was analysed in [14] where it was shown that its study can be 
reduced to that of a regular Lagrangian L = \L 2 . This Lagrangian L is quadratic 
in velocities and the dynamical vector field IY solution of the dynamical equation 
i(IY)cuL = <£El = cCL is not only a second order differential equation vector field 
but, moreover, it is a spray [22], the projection onto M 3 of its integral curves being 
the geodesies of the Levi-Civita connection defined by g. In other words, IY is 
the geodesic spray given by 

where the Christoffel symbols P jfc are 



■pi 1 „il 

1 jk ~ -j9 



dgki dgji dg jk 



dx J dx k dx l 
with g lJ being the inverse matrix of gij. 

In the particular case we are considering where g(v, w) = n 2 go(v, w), and using 
cartesian coordinates, 



P - 1 



dn dn „■ dn 
-51 + 



6< 



jk 



(5.3) 



dxi~ K ' dx k 3 dx 1 
It was also shown in [14] that the kernel of ujl is two-dimensional and it is 
generated by r L and the Liouville vector field A. The distribution V defined by 
Ker lvl is integrable because ujl is closed; actually [A, IY] = I\ and the distribu- 
tion is also generated by A and K defined by K = ^r L , for which [A, K] = 0. In 
cartesian coordinates K is expressed as follows: 



_d_ 

dx 1 



2 , . „ . \\v\\ 2 dn\ d 
-v l (v.Vn) - !LJ ^ 7 r^ 7T7 
n n ox 1 1 av l 



(5.4) 
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The theory of distributions suggests us the introduction of new local coordi- 
nates y % = F l (x, v), % = 1, . . . , 6, adapted to the distribution V defined by Ker ul, 
i.e., such that K = A = The search for these new coordinates is based 
on the solution of the partial differential equation system 

KF 1 = 1, AF 1 = 0, KF 2 = 0, AF 2 = 1, 

and 

RF 2+a = Qj AF 2+a = Q, for O = 1 . . . , 4. 

The explicit computation of these functions depends very much on the choice of 
the function n(x 1 , x 2 , x s ). We will illustrate next the theory with some particular 
examples. 

If n only depends on x 3 , the dynamical vector field is 

„ ,- d 2 i o dn d 2 9 o dn d 
T L = v 1 —- v V— r— v z V i 



dx 1 n dx 3 dv 1 n dx 3 dv 2 ^ ^ 

1 , x 2 ^ 2 2 v 3 2 \ dn d 
n dx 3 dv 3 ' 

and a function F 1 solution of AF 1 = should be F 1 = f 1 (x 1 , x 2 , x 3 , u 1 , u 2 ) with 

1 2 

u 1 = % y u 2 = \, and then the condition KF 1 = 1 reads 

KF 1 = u 1 ^- + u 2 ^- + 211-1^ — ^- -u 2 — ^- 
dx 1 dx 2 dx 3 n dx 3 du 1 n dx 3 du 2 

_ I(V 2 _|_ u 2 2 _ + u 2#/\ = x 

n cfa; 3 d-u 1 <9-u 2 

and the adjoint system is 

dx 1 dx 2 dx 3 ndu 1 ndu 2 df 1 

V ~ ^ ~ T" ~ u^u 12 + u 22 + 1)5 ~~ ~« 2 (w l2 + w 22 + l)^ ~" 
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^From here we see that f 1 is f 1 = x 3 + (p(x 1 , x 2 , u l ,u 2 ) where ip is an arbitrary 
function of the first integrals of the adjoint system 



dx 1 dx 2 dx 3 



n du 



n du 1 



u 1 



1 u\u 1 ' + u 2Z + 1) u 2 (u lZ + u 2Z + l)- dx . 

We see that C\ = ^ is one of such first integrals and then C 2 = x 1 — ^x 2 is 
a second one. The two last first integrals are to be obtained from 



dx 2 dx 3 



n du s 



1 W 2 ((1 + C 2)^ + 1) _ 

and the solution depends on the concret choice for the function n(x 3 ). More 



specifically, we find 



C3 = n exp 



' r 2 

J C((i + Ci 2 )C 2 + i 



where we must replace after doing the computation the constant C\ by the quotient 
C\ = ^7. The computation of the integral is quite easy and leads to 



C 3 = n\ 



(1 + C 2 )u 2 
(1 + C 2 )u 22 + 1 



= n\ 



l 2 1 2 2 

u l2 + u 22 + l 



For every value of C3 we can express u 2 as a function of n, 



u 2 = 



C 3 



vV-csxi + c?)' 



and then the fourth first integral is given by 



C 4 = x' 



Jo 



^(n 2 (()-C 2 )(l + C 2 ) 
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In a similar way we can look for the function F 2 amongst the solutions of 
KF 2 = and AF 2 = 1. Solutions of KF 2 = are the constants of motion and 
therefore z\ = n 2 ^ 1 , z 2 = n 2 v 2 which correspond to the momenta in the two first 
axes directions, where the medium index is constant, will give first integrals for K. 
Same for zs = x\ — ^%2, corresponding to the quotient of the third component of 
the angular momentum (in Optics called skewness function or Petzval invariant) 
by the second component of linear momentum and finally Z4 = n 
which is but the energy of the Lagrangian L. It is quite easy to see that F 2 = 
log l-y/z^l is a solution of the system KF 2 = and AF 2 = 1. 

According to this, we will do the following choice for the new coordinates: 



x ~ x , 



V 




x 



3 



(5.6) 
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namely, the inverse change is given by 



c. 



2 2 , 

x = y + 



x 3 = y\ 



v/(n J (0-(1)(l + C?) 
C 3 
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y 



and if we recall the expression of ujl written in the form 



(5.7) 



u>l = dx 1 A d 



nv 



\/ v 1 " 2 + v 2 " 2 + V s2 



+ dx 2 A d 



nv 



\/ v 1 " 2 + v 2,2 + v^ 2 



+ 



l 2 , 2 2 , 32N 3/2 
V + V + V^ 



r ' (/f ^ A cte 3 + v 2 d ) A f/.r :s 



we will get the new coordinate expression 

3 



U L = 



+ 



A V 



dy 2 + d 



o >/(« 2 (C)-C^)(l + C?) 



Ac% E 



1 



1 / 2 4 2 5 2^ ^ 4 ^ 4 A d y 3 + y 5 d y 5 A d v 3 ] , 

and therefore, when developing these expressions we will find that some terms 
cancel and only remain 



uj l = d [ y 1 + ^ y 2 ) A dy 4 + dy 2 A dy\ 



(5.8) 
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which shows that 

e =/ + 4 2 = * - £ f „ ^ a ^ 



v 2 Jo V(n 2 (0-C'i)(l + C'f) 

£ 2 = x 2 - f X ° 3 *r 

o V(« 2 (0-C7|)(l + C?) 



(5.9) 



and the corresponding 

1 4 77/1,1 2 <S n7j2 /r- 

v 1= y* = -7= 2 2 2 , ^ = zv 5 = -^= 2 2 2 (5.io) 

y v 1 + v 2 + v 3 y v l + v 2 + v 3 

are Darboux coordinates for the symplectic form induced in the quotient space. 

Let us now consider the particular but important case case in which the refrac- 
tive index becomes constant out of a region. If for x 3 > L, the index n is constant, 
the above mentioned Darboux coordinates and £ 2 are 



« 2 V(n 2 -C|)(l + Cf)' V(« 2 -c 2 )(i+c9' 

up to a constant, and from 

C =- ° 3 =- 

1 v 2 ' y/( n *-cri)(i + crt) v 3 ' 

we see that the Darboux coordinates become 

1 vl 3 2 1,2 3 n?j2 /r \ 

x -^i , x -^x, ===, ===, (5.11J 

in full agreement with (4.9). Therefore, for an optical system such that the refrac- 
tive index depends only on x 3 and, furthermore, the region in which the index is 
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not constant is bounded, we can choose Darboux coordinates by fixing a x 3 out- 
side this region and taking Darboux coordinates for the corresponding problem 
of constant index. This justify the choice of coordinates as usually done for the 
ingoing and outgoing light rays in the corresponding constant index media, i.e. it 
shows the convenience of using flat screens in far enough regions on the left and 
right respectively, and then this change of Darboux coordinates seems to be, from 
an active viewpoint, a canonical transformation. 

Let us now consider a different and interesting particular case in which there 
exists a symmetry axis and n depends on the distance to this axis alone. In this 
case the ray paths describe perfect optical instruments [17]. We will use cylindrical 
coordinates (r,9,z) and then the metric has nonzero elements g rr = n 2 , 
r 2 n 2 , and g zz = n 2 . The nonzero Christoffel symbols are 



rr 



V r — 

1 ee — 



1 dn 

n dr 

1 d / 2 2\ r 2 dn 

o~r{ r n ) = —t 

2n z dr n dr 



r r = _i^; 

zz n dr 



r e 



r e 

1 9r 



1 1 dn 

r n dr 
1 dn 



-pz -pz 

rz ~ zr ~ n dr ' 



(5.12) 



and therefore IY is given by 



+ 



d d 


d 


r dr Ve 00 


oz 


2 dn 


2 v r ve' 


-r v v r ' 

n dr 


r 



rvg 



_d_ 

8v r 



d 



+ 



2 dn 

■-—v e v r 
n dr 



_d_ 

dv. 



(5.13) 



30 J.F. Carinena and J. Nasarre 

Here we will take K as being given by The search for coordinates adapted 

to the distribution Keru;^ are found as in the previous example. A particular 
choice is: 

1 2 V 9 2 V V r + r2v d + V z 3 

v z nv z 

and 



4 /. f kdr 5 /" rfr 6 i 2 

y =6- =, y = z - I y = logn t> 

JO r 2 / c 2 n 2_fci_ 1 Jo / c 2 n 2_fc£_ 1 



where c and are the constants of motion 

A; = y 1 =r 2 — and c = y 2 = -^—^ 6 - - 

The presymplectic form lul, 

I TIT' 

<jj L = d6 A d ( ^^==^== ) +dr Ad 

+ dz Ad 




nv r 



a/v 2 + r 2 fg + 



^,2 



a/v 2 + r 2 Wg + u| 

when written using the new coordinates will become after cancellation of some 
terms 

to L = dy 4 A d ( - nr \ =] ■ dy 5 A d ( - = I . (5. 1 I) 

V V V r + r2v 9 + V z J V V V r + r2v 6 + V l 



from which we see that Darboux coordinates for the reduced symplectic form in 
the quotient space are 

dr 



e = y' = e- f . kdr a A 2 = v" = z- f - f 

JO r 2^ / / c 2 n 2_^_ 1 Jo ^ 



(5.15) 



c 2 n 2 - £r - 1 
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and the corresponding ones 






nv z 



(5.16) 



^Jv 2 + r 2 v' 2 + v 2 ' 



^Jv 2 + r 2 v 2 + v 2 ' 



which reduces to (4.11) when n is constant. 



6. Outlook 



The theory here developed suggests the interest of the study of what happens 
for anisotropic media, because of the recent ineterest in the use of anisotropic 
optical material. The theory can be reexamined along similar lines, because the 
basic principle of the theory is still the celebrated Fermat's principle of least time 
(or optical time if reflection is also allowed). The only difference is that when the 
medium is not isotropic the refractive index of the medium (given by the quotient 
n = |), may depend on the direction of the ray, and then An = 0. 

Again we will use the analogous mechanical problem where the Lagrangian 
function is given by L = n y/g(v, v). This Lagrangian is still homogeneus of 
degree one and the corresponding energy vanish identically, and therefore the 
Lagrangian is singular. It is possible to show that, at least when the Lagrangian 
L = |n 2 g(v, v) is regular, the curves solution for the Lagrangian L are just the 
curve solution of the original problem, even if the curves are reparametrized. Even 
more, the space of light rays in this particular case is also a symplectic manifold. 
This problem and the search for appropriate Darboux coordinates will be examined 
in a forthcoming paper. 



32 J.F. Carinena and J. Nasarre 

Acknowledgements: Partial financial support from DGICYT, Research Grants 
PS93-0582, is acknowledged. 



References 

[1] E.J. Saletan and A.H. Cromer, Theoretical Mechanics (J. Wiley, 1971) 

[2] R. Abraham and J.E. Marsden, Foundations of Mechanics (Reading, Ma: Ben- 
jamin, 1978) 

[3] V. Guillemin and S. Sternberg, Symplectic Techniques in Physics (Cambridge 
U.P., London, 1984) 

[4] P. Bamberg and S. Sternberg, A course on Mathematics for Students of Physics 
(Cambridge U.P., London, 1988) 

[5] J. Sanchez Mondragon and K.B. Wolf, Lie Methods in Optics , Lect. Notes in 
Physics 250, (Springer, 1986) 

[6] K.B. Wolf Lie Methods in Optics II , Lect. Notes in Physics 352, (Springer, 
1989) 

[7] A.J. Dragt and J.M. Finn, "Lie series and invariant functions for analytical 
symplectic maps", J. Math. Phys. 17, 2215-27 (1976) 

[8] A.J. Dragt, "Lie algebraic theory of geometrical optics and optical aberrations" , 
J. Opt. Soc. Amer. 72, 372-79 (1982) 



Symplectic structures in Optics 33 

[9] A.J. Dragt, "Lectures on nonlinear Orbit Dynamics" in Physics of High Energy 
Particle Accelerators , A. I. P. Conferences Proceedings No. 87. (Amer. Inst. 
Phys., New York, 1982) 

[10] A.J. Dragt and E. Forest, "Computation of nonlinear behavior of Hamiltonian 
systems using Lie algebraic methods", J. Math. Phys. 24, 2734-44 (1983) 

[11] K.B. Wolf, "Symmetry in Lie optics", Ann. Phys. 172 1-25 (1986) 

[12] O. Castanos, E. Lopez Moreno and K.B. Wolf, "Canonical transforms for 
paraxial wave optics", in Lie Methods in Optics , J. Sanchez Mondragon and 
K.B. Wolf, Lect. Notes in Physics 250 (Springer, 1986) 

[13] A.J. Dragt, E. Forest and K.B. Wolf, "Foundations of a Lie algebraic theory 
of geometrical optics" , in Lie Methods in Optics , J. Sanchez Mondragon and 
K.B. Wolf, Lect. Notes in Physics 250 (Springer, 1986) 

[14] J.F. Carinena and C. Lopez, "Symplectic structure on the set of geodesies of 
a Riemannian manifold", Int. J. Mod. Phys. 6, 431-44 (1991) 

[15] J.F. Carinena, "Theory of Singular Lagrangians" , Forts, der Phys. 38, 641- 
680 (1990) 

[16] J.E. Marsden and A. Weinstein, "Reduction of symplectic manifolds with sym- 
metry", Rep. Math. Phys. 5, 121-30 (1974) 

[17] D.D. Holm and G. Kovacic, "Homoclinic Chaos for Ray Optics in a fiber", 
Physica D 51, 177-188 (1991) 



34 J.F. Carinena and J. Nasarre 

[18] D.D. Holm and K.B. Wolf, "Lie Poisson description of Hamiltonian ray optics" , 
Physica D 51, 189-199 (1991) 

[19] D.D. Holm and G. Kovacic, "Homoclinic Chaos from Ray Optics in a fiber: 
200 Years after Lagrange", in Mechanics, Analysis and Geometry: 200 Years 
after Lagrange , Ed: M. Francaviglia (Elsevier, 1991) 

[20] K.B. Wolf, "Nonlinearity in Aberration Optic", in Symmetries and Nonlinear 
Phenomena , Eds: D. Levi D. and P. Winternitz, (World Sci. Press, 1988) 

[21] J.F. Carinena, J. Fernandez-Nunez and E. Martinez, "A geometric approach 
to Noether's Second Theorem in time-dependent Lagrangian Mechanics" , Lett. 
Math. Phys. 23, 51-63 (1991) 

[22] M. Crampin M. and F.A.E. Pirani, Applicable differential geometry (Cam- 
bridge U.P., London, 1986) 



